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Abstract. Rclativistic thermodynamics is derived from kinetic equilibrium 
in a general frame. Based on a novel interpretation of Lagrange multipliers 
in the equilibrium state we obtain a generic stable but first order relativistic 
dissipativc hydrodynamics. Although this was believed to be impossible, we 
circumvent this difficulty by a specific handling of the heat flow. 



1. Introduction 

Thermodynamics, hydrodynamics and kinetic theory are intimately related to 
each other. One may derive hydrodynamics from kinetic theory (via the Chapmann- 
Enskog expansion, or by momentum methods), but a certain interpretation of the 
thermodynamic quantities and relations is required prior the derivation. One may 
obtain the Gibbs relation by a proper averaging in the framework of a dynamic 
theory for continuous media. At the same time, a definite concept of local equilib- 
rium is necessary to calculate the entropy production and by doing to separate the 
dissipativc and non-dissipativc parts in the description. In this paper we explore 
these mutual relationships while generalizing rclativistic thermodynamics. By a 
novel treatment of Lagrange multiplicators, defining the local equilibrium state, we 
derive a first order and still generic stable, relativistic dissipativc hydrodynamics. 

Relativistic dissipativc hydrodynamics is not a straightforward generalization of 
the corresponding nonrelativistic theory: the simplest relativistic analogue of the 
Fourier-Navier-Stokes system of equations pQ is to date considered as conceptually 
wrong due to two reasons: it is unstable and acausal [2]. Thereafter generalized 
theories has been constructed by different methods, like gradient expansion in the 
Landau-Lifshitz frame H [5], GENERIC [6j [7], with the help of kinetic theory 
[HI El [101 [TH HSl [14], by phenomenological extensions [HI US E] or directly 
constructing divergence type equations |18j . The best known generalization, the 
Israel-Stewart theory [151 E El QH , extends the local equilibrium of Eckart theory 
improving its stability and causality properties [TH [20] . In spite of the motivation 
of the extension, there is little known about the stability and causality of other 
extended theories with the famous exception of the divergence type theories [2TJ [22j 
123] . These are namely causal by construction. In all the previously mentioned cases, 
the the Fourier-Navier-Stokes system is the base, which is extended relativistically 
by the Eckart theory, and extended further by the different methods listed above. 

From the thermodynamic point of view, the stability of the thermodynamic equi- 
librium of the fluid is indispensable: in the absence of external forces dissipation 
always damps perturbations. This property is called generic stability. However, 
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according to the investigations of Hiscock and Lindblom [21 [24] the Eckart theory 
is unstable in this sense, while the Israel- Stewart theory is conditionally stable. 
The appropriate conditions are far from being pure thermodynamical, their math- 
ematical form is complicated and lack a simple physical interpretation. The non- 
rclativistic Fourier-Navicr-Stokcs equations are generic stable, and the conditions 
of stability are straightforward: concave entropy function and the nonncgativity 
of the transport coefficients. In nonrclativistic hydrodynamics the first order the- 
ory is generic stable, while no first order generic stable theory was known in the 
relativistic case. This diametric contradiction in the stability properties between 
nonrclativistic and relativistic fluids and the complexity of stability conditions of 
the extended theories indicate that the thermodynamic background is insufficiently 
understood. 

The apparent departure between nonrclativistic and relativistic theories initi- 
ated several investigations regarding the possible reasons, not only to extend, but 
also to improve the Eckart theory. Garci'a-Cohn and Sandoval- Villalbazo suggested 
to introduce a balance equation of internal energy in addition to the balance of 
energy- momentum [25 [ I26 j . Muschik and Borzeszkowski noted, that this sugges- 
tion is not compatible with the foundations of relativity. We share their view, the 
energy- momentum density tensor has to include all energies of the system [2"7I |2"B] . 
Tsumura, Kunihiro and Onoshi [2H [301 131] used the renormalization group method 
in order to reduce the relativistic Boltzmann equation for describing slow macro- 
scopic dynamics. They concluded with the generic stability of their first order 
relativistic dissipative hydrodynamics in case of rarefied ideal gases in the Eckart 
frame. The obtained set of of dissipative hydrodynamic equations is similar for 
Garci'a-Cohn et. al. and for Tsumura et. al. in the sense that heat flow is propor- 
tional to the gradient of the temperature like in the nonrclativistic Fourier's law, 
the acceleration related term of the Eckart theory is missing. In our approach this 
infamous term is not left out, but other terms compensate their destabilizing ef- 
fect. The conditions for stability by missing accelerating term require the inequality 
dejd[i\ T > 0, which may be valid for uniform ideal gases, but can be violated in 
general, e.g. for Van der Waals gases. 

2. Eckart theory 

Here and in the followings wc use the metric diag(l, —1,-1,-1), the units are 
chosen such that the speed of light and the Boltzmann constant are unity c = 1, 
fee = 1. In Eckart theory the thermodynamic description of the local equilibrium 
is simple and reasonable: one introduces the entropy density as a function of the 
local rest frame energy density e = u a T ah Ub and of the particle number density 
n = u a N a . We note that in relativistic particle systems one of several conserved 
currents replace the particle number four- vector 7V a . The traditional continuity 
equation rather expresses the conservation of energy to leading order for massive 
(atomic) matter, eo = mc 2 n being the rest frame energy. Thereafter u a denotes the 
four velocity field actually used and T ab is the energy-momentum density tensor. 
The Gibbs relation in the Eckart theory is given as 

de = Tds + fidn, (1) 

where T is the temperature, s — u a S a is the entropy density, (i is the chemical 
potential and S a denotes the entropy four-current. 
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Analyzing the structure of the second law inequality, we recognized that the 
entropy must be at least a function of the energy-momentum density vector E a = 
T ab u b , instead of the energy density e [3.3 [33]. The particular case when the 
entropy depends on the Lorentzian length of the energy-momentum density vector, 
y/E a E a = \J e 2 + q a q a leads to the following form of the Gibbs relation: 

de + —dq a = Tds + fjdn, (2) 
e 

where q a = u c T cb A ba is the momentum density, the time-spacelike part of the 
energy-momentum density tensor with respect to the local velocity field, and A ba = 
5 ba - u a u b . 

The analysis of the structure of generic instabilities is also instructive in this 
approach. The energy-momentum density tensor of a single component fluid in a 
rest frame can be written as 

b= (u a T ab u b u c T cb A«\ f e q b \ 

\A^T cb u b A a c T cd A b d ) \q a P ab ) ' 1 > 

where the space-spacclike part P ab is the pressure, the current density of the mo- 
mentum. For one component, neutral fluids with zero spin, the energy-momentum 
tensor is symmetric. As a consequence, the time-spacclikc and space-timclike com- 
ponents are equafj]. However, the physical role of these spacelike vectors, q a and 
q b , is different. One of them is a momentum density, the other is an energy current 
density in a given frame. In nonrelativistic theories only the second one is related 
to dissipation, the pure momentum density is not. One may observe, that the in- 
stable modes of the Eckart theory are related to the momentum density. It is well 
known that the nonrelativistic approximation breaks the symmetry of these terms, 
momentum density and energy current density have different magnitudes when the 
speeds are low compared to the speed of light. However, the mentioned instabilities 
indicate a different hierarchy between the relativistically identical terms. The par- 
ticular form of the local equilibrium expressed by ([3 1 breaks the symmetrical role of 
the momentum density and energy current density in dissipativc rclativistic fluids. 
We note that Eq. ([3 1 already improves the stability properties of the corresponding 
a relativistic dissipative fluid mechanics, the theory proves to be generic stable not 
only in the Eckart frame but also in case of general flows. Moreover, the generic 
stability does not require conditions beyond natural thermodynamic requirements: 
the thermodynamic stability and nonncgative transport coefficients J33j [34] . 

A third, independent justification for a generalized Gibbs relation like J2J above 
arises investigating pure thermodynamic aspects. Analyzing the thermal interac- 
tion among two thermodynamic bodies, one gains a physical interpretation of the 
different transformation properties of the temperature suggested and discussed by 
Planck and Einstein, Blanusa and Ott, Landsberg [35] [33 [371 [33 [33 gU] , later in- 
corporating also the Doppler effect [JT] 03 H3 HI] • Since in our approach the Gibbs 
relation is modified at the density level, our explanation is more natural than those 
propositions where the transformation properties would depend on the definition 
of the relativistic thermodynamic body [45] H3 EH [48] . 

So far second law compatibility, generic stability of the fluid theory and the 
reasonable thermodynamic structure of the local equilibrium are only necessary 
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conditions to obtain a proper fluid theory. A further aspect arises from the com- 
patibility with the kinetic theory of gases. In the following we explore from the 
point of view of kinetic theory, whether the generalized rclativistic Gibbs relation 
([2]) can be further improved. Elaborating such an improvement we derive the corre- 
sponding rclativistic dissipativc hydrodynamic theory, which preserves all the good 
stability properties of the previous generalization always in the first order approach. 



3. Thermodynamic equilibrium in kinetic theory 

For the sake of simplicity, we consider in this paper a single-component Boltz- 
mann gas without external fields. The mass of the particles is m, their four- 
momenta are denoted by p a , \Jp a p a = tti. The Boltzmann equation for the one 
particle distribution function f(x a ,p a ) is given in the classical form 

P a d a f = C(f). (4) 

Here d a denotes the partial derivative of / regarding the spacetime variable x a 
and C(f) is the following collision integral: 

C(x, P ) = i J dwxdt/dwj [/7W,pib,Pi) - ffiW( P , Pl \p', P [)] , (5) 

where dui = is the momentum space measure. W(p,pi\p' ,p[) denotes the tran- 
sition rate from a particle pair momentum (p',p'i) before the collision, to particle 
pair momentum (p,pi) after the collision [49] . 

The particle number density four-vector, the energy-momentum density tensor 
and the entropy density four-vector arc defined as 



N a := J dujp a f, (6) 
T ab := J dujp a p b f, (7) 

S a := - f dup a f(]nf-l). (8) 

The local thermodynamic equilibrium is characterized by the condition of con- 
served entropy ([49], page 43). Then it follows, that the logarithm of the equilibrium 
distribution function is a collision invariant, therefore 

fo(x a , P b ) = e Q -^ c (9) 

with arbitrary fields a(x a ) and /3 c (x a ). 

For the equilibrium distribution function the particle number balance gives: 

daK = J dujp a d a f Q = Kd a a - T Qb a Q /3 b = 0. (10) 

In order to obtain thermodynamic relations, one can calculate the equilibrium value 
of the Boltzmann-Gibbs entropy (J5J): 

SS = J dup a f a (ln fo - 1) = (1 - a)N$ + p b T« b . (11) 

([TO]) and (33) result in 

d a SS + ad a NS - p b d a TZ» = 0. (12) 
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The formulas (|10[) and (|T2|) resemble a Legendre transformation, defining this way 
the thermodynamic structure of the theory with the five Lagrange multipliers a 
and Pb ISO]. However, (fT2)) cannot be interpreted as a generalized Gibbs relation. 
This relation alone docs not fix the entropy as a function of energy-momentum and 
particle number density. To find thermodynamic relations, one introduces a flow, 
a velocity field over the continuum material. 

According to the usual interpretation of the equilibrium distribution function 
©, such a vector field appears through the vectorial invariant (3 a , choosing that 
the velocity field parallel to it. However, in the general case, /3 a decomposes into 
components parallel and perpendicular to the velocity field 

r-^. (i3> 

There we have used the notation 1/T := f3 a u a and w a := A^/3 b /f3 c u c . It is 
convenient to introduce g a = u a + w a as well. We can see that g a u a = 1 and 
w a u a = 0, moreover, w a w a < 1, because f3 a , and therefore g a , are timelike. With 
/i = aT, the equilibrium distribution can be written as 

fo = e * . (14) 

The equilibrium particle number density four- vector and energy-momentum den- 
sity tensor can be calculated recognizing that the usual velocity field, defined by 
the direction of (3 a , is 



where w 2 — ~w a w a - In the following this usual equilibrium velocity field ii a is 
referred to as the natural frame. 

In the general frame, the equilibrium particle number four- vector TVq and energy- 
momentum density tensor Tq b are given as 

iV a = ng a = nu a + nw a , (15) 
TS b = (e+p)g a g b - P 6 ab 

n a„b 

a b , a b , b a \ ab , " " /i /?"\ 

= euu+qu+qu — pA H . (16) 

e + p 

Here, the thermodynamic quantities in the general frame (without hat) are related 
to the usual quantities in the natural frame (with hat) as T = Vl — w 2 T, fj, = 
\/l — w 2 fi, n = h/ \/l — w 2 . p = p, e= (e + pw 2 )/(l — w 2 ), (see [51| for more 
details). The relation between the momentum density and the particle number 
current density is fixed: 

q a = (e+p)w a . (17) 

Therefore, for nondissipative fluids, the Eckart frame and the Landau-Lifshitz 
frame coincide. Substituting (fTS"!) and (jTB")) into (jTTJ) yields 

SS = ^{nT- f xn + E b g b ), (18) 

where E a = eu a + q a . The entropy current density in the general frame is paral- 
lel to g a . The thermodynamic relations are calculated by requiring the following 
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inequality (5(3 EQ: substituting (jUJ and into (JT2 



9 a 5 a + ad a N« - Ma?? = 

9 — (Td a s + ad a n - g b d a E b - pw b d a u b ) + ( Ts + M n - . 96 £ fa -p)=0. (19) 

Now the first term of the right hand side of the previous equality does not 
contain contracted derivatives, and allows to extract thermodynamic relations. One 
recognizes that the entropy density in the general frame is the function of the 
particle number density, of the energy-momentum density four-vector and of the 
velocity s = s(E a , n, u a ), with the following partial derivatives: 

Hi = =H ds = il — = OHl (90) 

On a T" dE a T ' du a T 1 ' 

The consequent Gibbs relation for the differentials is 

g a dE a + pw a du a = de + w a dq a + [(e + p)w a - q a ]du a = Tds + ^dn. (21) 

Utilizing Eq. (|T7l) . the Gibbs relation defining the local equilibrium of a rela- 
tivistic fluid is given by 

de + j-dq a = Tds + /j,dn, (22) 

where h = e + p is the enthalpy density. Whenever w a is zero, we obtain JIJ as a 
special case. This unusual Gibbs-relation already considers both energy-momentum 
exchange in the form of transport and also in the form of mechanical work. The 
Legendre transformation properties are fixed by the last bracket in (|19|) . It vanishes 
if 

p + g a E a =p + e+ ^ =Ts + fin. (23) 
h 

A comparison to (fT5| provides the thermal equation of state of our ideal gas: 
p = nT. This equation of state is valid in general frames. 



4. DlSSIPATIVE HYDRODYNAMICS 

In order to determine the constitutive functions of the hydrodynamic theory in a 
linear Onsagcrian framework, we calculate the entropy production with the help of 
the generalized Gibbs relation (f2"2"j) . Introducing the substantial or comoving time 
derivative ^ := u a d a , denoted by an overdot, the basic balance equations are 

d a N a - n + nd a u a + d a f = 0, (24) 
d b T ab = eu a + eu a d b u b + q a + q a d b u b + eii a 

+ u a d b q b + q b d b u a + d b P ab = a . (25) 

The energy and three-momentum balances are the time- and spacelike compo- 
nents of the energy-momentum balance: 

u a d b T ab = e + ed b u b + u a q a + d b q b + u a d b P ab = 0, (26) 
A a c d b T cb = eu a + A^q b + q a d b u b + q b d b u a + A a c d b P cb = a . (27) 

The calculation of the entropy production is now straightforward. The time de- 
rivative of the entropy density is determined by using the functional dependencies 
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inherent in the Gibbs relation (f22[) . Then we substitute the balance equations (|24|) 
and (|26j) and we obtain 

d a S a = s + sd b u b + d a J a + d a ( ^ a ~ T ^ a 
= ± (P ab ~{e + q -f-Ts- ,n) A-) 8 a u b + <f * + | + 

= -fdaj; + ^Tl ab d a u b + q a X? h > 0. (28) 

Here the static pressure was determined by ([2"B"| , LT ab = P ab + pA ab is the viscous 
stress and 

Xi = A ba (d a ^ + -^ (hu a + q a + q a d c u c )\ (29) 

is the thermodynamic force of the thermal interaction. The classical form of the 
entropy current density J a = ^(q a — nj a ) in ([2^)1 can be justified by a deeper 
analysis of the second law [3"2"| . 

Consequently assuming linear relationship between the thermodynamic fluxes 
and forces leads to the following constitutive relations: 

q a = \X? h , f = -iA ab d b ^ IF 6 = -2r t d (a u u > - c u c A ab , (30) 

where A is the heat conduction coefficient, £ is the diffusion coefficient, r), C are 
the shear and bulk viscosities, respectively. The bracket denotes the symmetric 
traceless part of the spacelike tensor 

d (a u b) = A a A b ( d C U d + d d U C _ l g d 



The remarkable part of the entropy production formula is (|29|) . The gradient of 
the inverse temperature and the acceleration solely constitute the driving force 
of thermal interaction in the Eckart theory. The additional terms containing the 
momentum density q a are absent there. The terms beyond the gradient of the 
temperature stem from the first part of the momentum balance (|27j) divided by the 
enthalpy density. This indicates that by relativistic heat conduction the momentum 
and energy transport cannot be separated. 



5. Generic stability 

A necessary, expected requirement of fluid theories is generic stability, that is 
the linear stability of the homogeneous equilibrium of the system of equations (|24p , 
(|2l)]) . (f2"T|) and (f5U|) . In the following we sketch the proof of this property in the 
Eckart frame, where j a = a . 
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Denoting the perturbed fields by S we can put the linearized equations (|24|) . (|26l) . 
27| and the first and last equations of ([30]) into the following form 

= Sn + nd a Su a , (31) 

= 5e + hd a 6u a + d a Sq a , (32) 

= hSu a + A ab d b 6 P + Sq a + A a c d b SU cb , (33) 

s 1b\ 



= 



= 




1 /ifaf, + 

Obdy y 



(34) 



i)A oc A M (9 c fo d + d <Ju c - -d f Su f A cd ). (35) 



Now we are looking for solutions in the form of exponential plane wave form SQ = 
Q e Tt+ ' lkx . Then the condition of asymptotic stability of the above equations is 
that the solutions of the characteristic equation detM — (dct N)(dct R) 2 = give 
negative real parts for T. Here 



(36) 
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Here h = e + p, fj - 


= ^ 


+ 1^7, »7 



(37) 



3'' 

fjk 2 



0. 



A = X/T 2 . The determinant of R gives the condition 
which results in negative real T solutions for all wave numbers k. The determinant 
of N results in the following condition 



T A \h + k 2 



, AT. 



-V 



2;„2 



T Z k 



I) 



Xnd n T 



XT 



nd n p 



Tk 2 [hd e p + nd n p + k 2 f)Xd e T) + k 4 Xn {d n pd e T - d e pd n T) = 0. 

According to Routh-Hurwitz criteria the real parts of the roots of the polinomial 
equation aox 3 + a\x 2 + a,2X + 0,3 =0 are non negative if the coefficients are all 
non-negative and a^a-i — a^a^ > [52]. It is straightforward to check that these 
inequalities arc fulfilled, if the transport coefficients A, 77, ^are nonnegative and the 
inequalities of the thermodynamic stability are satisfied, i.e. d e T > 0, <9„#r > 



and also d e Td n % ~ (%M 



> 0. 



One may prove similarly the generic stability of the above hydrodynamic equa- 
tions in case of general frame, too. 



6. Summary and conclusions 



In nonrelativistic mechanics three-momentum and three-velocity are parallel, 
in relativistic mechanics the corresponding respective four-vectors are orthogonal. 
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While in nonrelativistic mechanics the momentum is always zero in a local rest 
frame - from v = follows mv = , in a relativistic theory it is not necessarily 
so, there is only a particular rest frame where the momentum vanishes. Therefore, 
in case of thermal interactions between at least two participants, one of them will 
see a nonzero exchanged momentum. That is the basic reason why relativistic and 
nonrelativistic thermodynamics are different, and that is what a proper relativistic 
concept of local equilibrium should reflect. It is still a local equilibrium by definition, 
since the entropy production is zero, and the generalized Gibbs relation requires 
a constant velocity of the thermodynamic body. However, it depends on the rest 
frame momentum because the relativistic concept of local equilibrium should be 
formulated in a general frame. 

In this paper we have shown that a particular extension of the concept of local 
equilibrium is compatible with kinetic theory when considering general flows. A 
corresponding Gibbs relation and the related thermodynamic formalism introduced 
the rest frame momentum density q a as an additional density variable and the 
related intensive quantity w a = q a /h is a particular spacclikc quantity with the 
physical dimension of velocity. The necessity of the distinction between momentum 
density and energy current density in a dissipative theory was already recognized in 
[55] and Eq. ^ was suggested as the basic formula defining the local equilibrium. 
Our present study of relations to kinetic theory revealed that in local equilibrium 
inertia is defined through the balance of momentum, therefore, the enthalpy density 
substituts the energy density in the properly generalized Gibbs relation ([22]) . 

From the point of view of the traditional concept of local equilibrium by Eckart, 
our approach might be considered as second order in the heat flow. According to 
the definition given by Hiscock and Lindblom [2] the obtained relativistic dissipative 
hydrodynamics is first order, because the deviations from the local equilibrium are 
first order in the entropy current. 

The generic stability of the equations is proved in the Eckart frame. 
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